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Finite-temperature phase transitions in quasi-one-dimensional quarter-filled systems are in- 
vestigated using the extended Hubbard model with electron-lattice coupling. By a quantum 
Monte Carlo method combined with interchain mean-field approximation, we clarify compet- 
ing and coexisting behaviors among charge ordering, dimer Mott, and spin-Peierls states. It 
is pointed out that an anharmonicity of lattice distortions plays an important role in multi- 
critical behaviors. The results are compared with experimental data for quasi-one-dimensional 
molecular conductors such as DCNQI and TMTTF compounds. 
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Quasi-one-dimensional (Q1D) molecular conductors 
have been recognized as materials suitable for studying 
the effects of strong correlation. In these compounds, 
Coulomb interaction and enhanced fluctuations due to 
low dimensionality give rise to keen competitions in 
charge, spin, and lattice degrees of freedoms. As a re- 
sult, they show various phase transitions despite the fact 
that their nonintcracting band structures are often very 
similar. 1 

A typical case is the family of (i?ii?2-DCNQI)2X (Ri, 
i?2 : substituents, X: monovalent cation). The electronic 
structure is commonly described by a Q1D quarter-filled 
7r-band of DCNQI chains, 2 whereas their physical proper- 
ties depend strongly on R\, R2, and X. For instance, DI- 
DCNQI 2 Ag shows a charge-ordering (CO) phase tran- 
sition at T = 220 K and an antifcrromagnetic transi- 
tion at 5 K. 3 ' 4 On the other hand, DMe-DCNQI 2 Ag ex- 
hibits, instead of CO, lattice dimerization at T = 100 
K, which drives the system effectively half-filled, result- 
ing in a dimer Mott (DM) insulating state. Furthermore, 
tctramcrization occurs at T — 80 K, which is ascribed to 
a spin-Peierls (SP) transition. 5 Such a variety of prop- 
erties suggests a subtle balance among different phases 
under electron correlation. 

Another good example is found in (TMTTF) 2 A (X: 
monovalent anion), whose Q1D 7r-band of TMTTF 
chains is quarter-filled in terms of holes. 6 In contrast 
to the DCNQI systems, an intrinsic dimerization exists 
along the chains from the outset. These compounds also 
exhibit rich phases depending on X. The CO transition 
occurs at around 100 K, which is followed by a SP tran- 
sition at a lower temperature (T) for X = PF@ and AsFg 
or an antifcrromagnetic transition for X= SbFg. 7 ~ 9 For 
X = AsF 6 , the CO phase is suppressed under applied 
pressure (P), while the SP state persists up to higher P. 
This indicates that the CO and SP states are competing 
in nature. 8 



Theoretically, correlation effects in such Q1D molecu- 
lar conductors have been studied using the ID or Q1D 
quarter-filled extended Hubbard model with on-site and 
intersite Coulomb interactions, mainly focusing on the 
ground state. 10 Several studies have recently been con- 
ducted to describe their finite-T properties. n ~ 16 In such 
studies, it is crucial to treat ID electronic fluctuations 
properly, which originate from the strong electron corre- 
lation. This is because development of correlation lengths 
is strongly anisotropic along ID chains, which finally re- 
sults in phase transitions through interchain interaction 
or coupling to the three-dimensional lattice. Moreover, 
the paramagnetic insulating nature in strongly corre- 
lated phases such as the CO and DM states can be cap- 
tured only when the fluctuation effects are taken into 
account. Recently, a numerical study including such ID 
fluctuation effects was performed 14 using the quantum 
transfer-matrix (QTM) method combined with the inter- 
chain mean-field approximation. Although competitions 
among different phases were elucidated, the low-T region 
was not fully investigated because of limitations of the 
QTM method. Most of the intriguing competitions in real 
materials appear at much lower T; hence, it is strongly 
desired to study the low-T region to understand the ex- 
perimental results. 

In this paper, we investigate the finite-T properties of 
Q1D quarter-filled systems, particularly finite-T phase 
diagrams including various charge and lattice ordered 
phases, by a quantum Monte Carlo method with inter- 
chain mean-field approximation. This numerical method 
is reliable enough to obtain high-precision results down 
to low T far below that previously studied by QTM. The 
main results are summarized in the phase diagrams in 
Fig. 1, which show how phase competition and coexis- 
tence are controlled by electron correlation and electron- 
lattice coupling. We will discuss the results in comparison 
with experiments on the above-mentioned compounds. 
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Fig. 1. Phase diagrams in the plane of zV± and T for t = 1, U = 6, V = 2.5, and Kp = 0.75; (a) 8^ ~ and 
Kp 2 = 0, (b) <5d = and Kp 2 = 0.75, (c) <5d = 0.02 and Kp 2 = 0. Phase boundaries at T = (cross symbols) 
are obtained by the exact diagonalization (Lanczos) method for iV=12. (d) Schematic views of different ordered 
states. Dashed ellipses represent spin-singlets. Note that in (c), intrinsic dimerization always exists; therefore 
CO(+iDM) and SP(+iDM) phases respectively have the same patterns as DM+CO and DM+SP states in (d). 
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Our Q1D Hamiltonian is given by H = £\. ^id +^ii 
where 7i iD and 7ix represent the intrachain Hamilto- 
nian of the j'-th chain and interchain one, respectively. 
The intrachain part is described by the extended Hub- 
bard model with the Peierls-type electron-lattice cou- 
pling, whose Hamiltonian reads 
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(1) 

where the notations are referred to ref. 14. The lattice 
distortions u,j couple to electrons through the transfer 
integrals, and the coupling constant is incorporated in 
the definition of Uij. We treat Uij as classical variables. 
The elastic energy is considered up to the fourth order, 
where the term with Kp 2 represents the anharmonic con- 
tribution. The interchain part is given by 



Hx = V±J~ y j E n ''J n '' fc ' 



(2) 



where V± denotes the interchain Coulomb interaction 
and the summation for (j, k) runs over nearest-neighbor 
chains. In the following calculations, we choose the on- 
site and nearest-neighbor Coulomb interactions as U = 6 
and V = 2.5, respectively, 17 and the elastic constant as 
Kp = 0.75, in energy unit of t. The chemical potential [i 
is controlled so that the system is at quarter filling. 

We derive an effective ID model under two approxi- 
mations. The first is the interchain mean-field approxi- 
mation; by extending eq. (4) in ref. 14 to incorporate the 
CO+SP state in Fig 1(d), we assume the mean- field form 
of H± as 

nT = ^ Yl Mm-!) + (n l+1 )) m - {m-iHn^i)} , 

i 

(3) 
where z denotes the number of nearest-neighbor chains 
and the chain index j is dropped hereafter. The second 
is the adiabatic approximation for the lattice distortions; 
we consider the situation where the dynamics of lattices 
is much slower than that of electrons and fluctuations of 
electrons play a major role. The lattice distortions Ui are 
then determined so as to minimize the free energy under 
the constraint Yli u i = 0- Note that we omit the chain in- 
dex j again, u,,j = Uj. These two approximations give rise 
to an effective ID model Tiip, + 7Y r f IF with self-consistent 
conditions for (rii) and itj, which exhibits finite-T phase 
transitions. This is suitable for our purpose, i.e., to cap- 
ture the essential features of the phase transitions gov- 
erned by ID electronic fluctuations which occur at finite 
T assisted by the additional three-dimensionality. 

We solve the effective ID model by a quantum 
Monte Carlo technique called the stochastic-series- 
cxpansion (SSE) method 18 ' 19 in the operator-loop- 
update scheme. 20,21 The SSE calculations fully include 
thermal and quantum fluctuations of electrons and give 



unbiased and high-precision data for the effective ID 
model. We use the expectation values by SSE in the 
self-consistent equations, and obtain the inputs for the 
next SSE calculations in turn. This cycle is repeated until 
the convergence is reached for all (rii) and Uj. We con- 
sider symmetry breaking with a four-site period along 
the chain by taking account of tetramerizations in the 
experiments, and solve the self-consistent equations for 
eight order parameters: four each for the charge densi- 
ties (ni) and the lattice distortions ui (I — 1,2,3,4). We 
have calculated systems with sizes up to N — 64 sites 
and confirmed that finite-size effects are negligible down 
to T = 0.02*. 

First, let us show the results in the absence of in- 
trinsic dimerization, i.e., <5d = 0. Figure 2 shows the In- 
dependences of the charge densities (rii), the lattice dis- 
tortions Ui , and the charge and magnetic susceptibilities 
Xc and Xsi hi the case of the harmonic lattice (Kp 2 = 0). 
Figures 2(a)-2(c) show the results for zV± = 1.5, which 
are typical behaviors when the interchain Coulomb in- 
teraction is dominant compared with the electron-lattice 
coupling, and Figs. 2(d)-2(f) show those for an opposite 
situation with zV± = 0.25. 

In the former case, CO takes place at T = Tqo — 
0.86, as observed in the alternation of (rii) shown in 
Fig. 2(a). The charge susceptibility Xc (= dn/d/j,, n: av- 
erage electron density) suddenly drops at Tco and de- 
creases rapidly below it owing to the opening of a charge 
gap, while the magnetic susceptibility Xs (= dm/dh, m: 
average magnetic moment, h: magnetic field) show no ob- 
vious change at Tco 14,22 [Fig. 2(c)]. This is the transition 
to a paramagnetic insulating state. The order parameter 
for CO develops down to T ~ 0.3, below which, unex- 
pectedly, it slightly decreases [the inset of Fig. 2(a)]. This 
temperature corresponds to the region where Xs devi- 
ates from the Curie- Weiss behavior [Fig. 2(c)], i.e., where 
spins start to interact with each other. Hence, we con- 
sider that this reduction of CO is due to ID fluctuations 
fully taken into account in our scheme. With further de- 
crease in T, tetramerization emerges in both the charge 
and bond sectors at T = Tco+sp — 0.06, as shown in 
Figs. 2(a) and 2(b), which was not accessible in the pre- 
vious QTM results. 14 This transition is ascribed to the 
SP transition, where two neighboring spins at charge- rich 
sites form a spin-singlet pair with spontaneous lattice dis- 
tortions: the low-T phase is the CO+SP state 17 sketched 
in Fig. 1(d). The sudden drop in Xs below T co +sp m 
Fig. 2(c) confirms that the system has a spin gap. 

Different successive transitions occur when the 
electron-lattice coupling is dominant [Figs. 2(d)-2(f)]. In 
this case, the lattice dimerization first takes place at 
T = Tdm — 0.67 [Fig. 2(e)], which leads the system effec- 
tively half-filled. This indicates the DM transition with 
the opening of a charge gap, as seen in the strong sup- 
pression of Xc [Fig. 2(f)]. Xs shows no significant change 
at TrjM) and again this is characteristic of the transition 
to a paramagnetic insulator. Similarly to that in the CO 
case, the DM order parameter decreases slightly below 
T ~ 0.2 [the inset of Fig. 2(e)], where x s deviates from 
the Curie- Weiss behavior [Fig. 2(f)]. Symmetry breaking 
with a four-site period occurs at T = Tdm+sp — 0.08, as 
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shown in Figs. 2(d) and 2(e). This is another SP tran- 
sition from the DM state to the DM+SP state 23 with 
a sudden drop in Xs [Fig. 2(f)], whose ordering pattern 
is different from that of the CO+SP state, as shown in 
Fig. 1(d). 
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ter of multicriticality. It is clear that the anharmonicity 
Kp 2 stabilizes the coexistence of DM and CO. We note 
that the harmonic model with Kp 2 = appears to be on 
the verge of bicritical and tetracritical behaviors, 24 even 
though it is rather difficult to exclude the coexistence 
in a very narrow range. Typical T dependences of the 
order parameters in the coexisting regime are shown in 
Fig. 3 for zV± = 0.9. At Tdm — 0.67, dimerization oc- 
curs first in the bond sector, and CO appears at a lower 
7dm+co — 0.42; CO coexists with DM below Tdm+co- 
With further decreasing T, the degree of lattice dimer- 
ization decreases as CO grows, and finally, the system 
exhibits a transition to the CO+SP state. This transi- 
tion appears to be of first order since the lattice dimer- 
ization suddenly vanishes as shown in Fig. 3(b). For large 
or small zV±, the coexistence of DM and CO disappears, 
and the successive transitions as in Fig. 2 take place. The 
results are summarized in Fig. 1(b). 
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Fig. 2. Temperature dependences of charge densities, lattice dis- 
tortions, and charge and magnetic susceptibilities for (7 = 6, 
V = 2.5, K P = 0.75, <5 d = 0, and Kp 2 = 0; (a)-(c) zV ± = 1.5, 
(d)-(f) zV± = 0.25. Open and closed circles represent results 
for JV=32 and 64, respectively. Statistical errors arc smaller 
than the symbol sizes. Data shown by cross symbols at T=0 
are obtained by the exact diagonalization (Lanczos) method for 
N=12. The dashed line in (c) [(f)] shows a Curie- Weiss fit for 
0.3 < T < T co (T DM ). 



Results obtained by varying zV± are summarized in 
the zV±-T phase diagram of Fig. 1(a). At high T, there 
are phase transitions from a high-T uniform metal to 
intermediate-T paramagnetic insulators with the two-site 
period, the DM or CO state; both transitions are of sec- 
ond order. The DM and CO states are competing with 
each other and the transition between them appears to 
be of first order. Thus, the phase diagram is likely of 
bicritical type as a result of this competition. We will 
come back to this point later. At lower T, the tetramer- 
izations emerge upon both the CO and DM states, and 
these SP transitions are of second order. Note that the 
low-T phases were not identified in the previous QTM 
study because of numerical limitations. 14 

The phase diagram is also investigated in the case of 
anharmonic lattice distortions, as shown in Fig. 1(b), for 
Kp 2 = 0.75. The anharmonicity opens a window where 
DM and CO coexist, and therefore, the phase diagram 
shows a tetracritical behavior. Although the coexisting 
phase was already found in ref. 14, here we elucidate 
that the lattice anharmonicity is a controlling parame- 



Fig. 3. Temperature dependences of (a) charge densities and (b) 
lattice distortions for U = 6, V = 2.5, K P = 0.75, <5 d = 0, 
Kp 2 = 0.75, and zV± = 0.9. The same symbols as those in 
Fig. 2 are used. 



Finally, we study the effects of the intrinsic lattice 
dimerization <5 d . The phase diagram for 5^ = 0.02 is 
represented in Fig. 1(c) and typical T dependences of 
physical quantities are shown in Fig. 4. For 5d ^ 0, there 
is no DM transition since the lattice dimerization is al- 
ways induced even at high T [Figs. 4(b) and 4(e)]. We 
thus refer to the dimerized state without spontaneous 
symmetry breaking at high T as the intrinsic DM (iDM) 
state in the phase diagram in Fig. 1(c). When zV± is 
dominant, CO and SP transitions occur similarly to the 
cases of 5d = as shown in Figs. 4(a) and 4(b). Note, 
however, that the CO+SP(+iDM) state here is slightly 
different from the CO+SP state for (5 d = because of 
the underlying lattice dimerization [Fig. 1(d)]. On the 
other hand, when the electron-lattice coupling is domi- 
nant, CO is suppressed and the iDM state prevails until 
the tctramcrization takes plate at lower T [Figs. 4(d) 
and 4(e)]. Xs behaves similarly to that for <5 d = in both 
cases [Figs. 4(c) and 4(f)], indicating that the intrinsic 
lattice dimerization have no significant effect on the spin 
sector. On the other hand, Xc behaves differently at high 
T: it is suppressed even at T > Tqo because of the Mott 
insulating nature of the iDM state. The phase diagram 
is summarized in Fig. 1(c), where the DM transition is 
absent and the phase transition between the iDM and 
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CO(+iDM) phases is of second order. 
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Fig. 4. Temperature dependences of charge densities, lattice dis- 
tortions, and charge and magnetic susceptibilities for [7 = 6, 
V = 2.5, K P = 0.75, 5 d = 0.02, and A'p 2 = 0; (a)-(c) zV ± = 1.2, 
(d)-(f) zV± = 0.2. The same symbols as those in Fig. 2 are used. 



Now let us discuss our results in comparison with the 
experimental data. The family of (i?ii?2-DCNQI)2A has 
uniform DCNQI chains at high T, corresponding to the 
cases of (5d = in our model. The successive transitions 
with dimerization and tetramerization observed in DMc- 
DCNQI 2 Ag 5 are reproduced by the DM and SP transi- 
tions, as shown in Figs. 2(d)-2(f). The CO transition in 
DI-DCNQI 2 Ag 4 is also reproduced in Figs. 2(a)-2(c), but 
the low-T antiferromagnetism 3 is out of the scope of our 
model because we neglect the interchain hopping that re- 
sults in an effective antifcrromagnctic exchange interac- 
tion between chains. 25 On the other hand, (TMTTF) 2 A 
has an intrinsic dimerization, and hence corresponds to 
Sd ^ 0. The CO and SP transitions observed for A"=PF 6 
and AsFg are reproduced in Figs. 4(a)-4(c). Furthermore, 
our phase diagram in Fig. 1(c) shows a good agreement 
with the P-T phase diagram for A=AsFg 8 when we read 
the interchain Coulomb interaction zV± as the inverse 
of P. This is reasonable since the pressure is supposed 
to increase mainly the intrachain transfer integrals, i.e., 
to effectively decrease the interchain interaction. As to 
the pressure effect on DCNQI compounds, the case for 
DMe-DCNQ^Ag has not been explored yet, but it was 
shown that DI-DCNQI 2 Ag exhibits a peculiar tricriti- 
cal behavior under P. 26 Although we have found no tri- 
critical point in Figs. 1(a) and 1(b), our observation of 
the role of lattice anharmonicity in the multicriticality 
indicates the possibility of capturing such behavior by 



extending our model, particularly in the electron-lattice 
part. 

It should be pointed out that the obtained phase di- 
agrams will be modified if we go beyond the interchain 
mean-field or adiabatic approximation. For example, the 
effects of fluctuations and quantum nature in the lat- 
tice are anticipated to contract the DM and SP phases. 
These phases would nevertheless survive at finite T pro- 
vided that three-dimensional lattice couplings are in- 
cluded. Hence, we presume that the general features of 
our phase diagrams remain robust, and bear the com- 
parisons with experiments. The effects of fluctuations 
neglected in the present calculations are left for future 
study. 

In summary, we have investigated the finite- 
temperature phase transitions in quasi-onc-dimcnsional 
quarter-filled systems using the extended Hubbard 
model with electron-lattice coupling. The effective one- 
dimensional model, obtained by interchain mean-field 
approximation, has been numerically solved by the 
stochastic-series-expansion Monte Carlo method down to 
temperatures far below that previously studied. A variety 
of phase diagrams have been explored and rich behav- 
iors of competition and coexistence have been clarified 
among charge, lattice and spin degrees of freedom. We 
have found that an anharmonicity of lattice distortions is 
a key parameter of multicritical behaviors. Our results re- 
produce well the charge ordering, dimer Mott, and spin- 
Peierls states in DCNQI and TMTTF compounds. 
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